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SUMMARY 


An analytical and experimental study was carried out to understand the exten- 
sional and flexural behavior of multiwall sandwich, a metallic insulation composed of 
alternate layers of flat and dimpled foil. The multiwall sandwich was structurally 
analyzed by using several simplifying assumptions combined with a finite element 
analysis. The simplifying assumptions made in this analysis were evaluated by bend- 
ing and tensile tests. Test results validated the assumption that flat sheets in 
compression do not significantly contribute to the flexural stiffness of multiwall 
sandwich for the multiwall geometry tested. However, calculations showed that 
thicker flat sheets may contribute significantly to bending stiffness and cannot be 
ignored. 

Both test data and results of a finite element analysis show that the dimpled 
sheet is stiffer in the 0° direction than in the 45° direction. Tensile tests show 
that in the 0° direction, the dimpled sheet behaves linearly until material yielding 
occurs; however, in the 45° direction the dimpled sheet exhibits nonlinear behavior 
as a result of significant shape changes. The nonlinear behavior is strongly depen- 
dent on the lateral-edge constraint of the dimpled sheet. Two sets of boundary con- 
ditions were used in the linear finite element analysis to bracket the actual con- 
straint of the dimpled sheet in the multiwall sandwich. Results of this analytical 
approach compare well with test data; both show that the extensional stiffness of the 
dimpled sheet in the 0° direction is about 30 percent of that for a flat sheet, and 
that in the 45° direction, it is about 10 percent. The analytical and experimental 
multiwall bending stiffnesses showed good agreement for the particular geometry 
tested. 


INTRODUCTION 

Multiwall sandwich is a unique structural sandwich concept, which was originally 
designed as a vacuum-sealed insulation for cryogenic tankage (ref. 1). As an out- 
growth from this concept, several studies have considered an unsealed version of the 
multiwall sandwich for use as a thermal protection system (TPS) (refs. 2 to 4) . An 
exploded view and cross section of a typical TPS tile made of multiwall sandwich are 
shown in figure 1 • The sandwich consists of alternate layers of flat and dimpled, 
foil-gauge metal sheets bonded together at the crests of the dimples. As shown in 
the cross-sectional view, bonding of the flat and dimpled sheets forms a complex, 
three-dimensional structure. 

Although designed primarily from thermal considerations to be a lightweight 
insulation, the multiwall sandwich must also be designed to carry the loads to which 
it will be subjected in TPS applications. Structural analysis of multiwall is com- 
plicated by buckling of the flat sheets loaded in compression and by the complex 
shape of the dimpled sheet. In the initial conceptual development of multiwall for 
TPS application, reference 3, approximate engineering-type calculations were 
developed for preliminary structural analysis. These calculations are based on the 
assumptions that (1) flat sheets in compression were totally ineffective, (2) dim- 
pled sheets were fully effective in both tension and compression, and (3) each 
dimpled sheet could be represented as a flat sheet located at the centroid of the 
dimpled sheet. 



The purpose of this paper is to study the structural behavior of multiwall sand- 
wich (referred to as multiwall) and to develop a more refined, but still approximate, 
method for calculating its flexural stiffness. Results from tensile tests of dimpled 
sheets and bending tests of multiwall-sandwich specimens are used to evaluate some of 
the analytical assumptions and to provide experimental data for comparison with 
analytical predictions. 
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SYMBOLS 

effective area per unit width of each sheet, used to calculate moment of 
inertia (see appendix B) , in 2 /in. 

distance from outer support to inner load application point in four-point 
bending tests (see appendix D) , in. 

arbitrary column width, used for buckling calculations (see appendix A), 
in. 

distance from neutral bending axis, in. 

material modulus of elasticity, psi 

effective modulus of elasticity, average stress divided by average strain, 
psi 

dimple height, the vertical distance between dimple peaks and troughs (see 
appendix B), in. 

bending moment of inertia per unit width, in 4 /in. 

length of the boundary on the finite element model, in. 

distance between inner load application points in four-point bending tests 
(see appendix D), in. 

bending moment, in-lb 

constant bending moment applied to center section of four-point bending 
specimens (see appendix D) , in-lb 

number of effective sheets in a multiwall sandwich under bending (see 
appendix B) 

applied load in bending tests (see appendix D) , lb 

ratio of the thickness of a flat sheet to that of a dimpled sheet with the 
same extensional stiffness 

distance between nodes along the 45° direction in dimpled sheets, in. 

dimpled-sheet thickness before forming, in. 

flat-sheet thickness, in. 
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face-sheet thickness, in. 
transverse displacement of beam, in. 

distance from lower surface of sandwich to neutral bending axis (see appen- 
dix B), in. 

measured displacement in bending tests, in. 
stress, psi 
Subscripts: 

ab after flat sheets buckle 

av average 

bb before flat sheets buckle 

cr critical 

0° flat sheet supported by a dimpled sheet oriented in the 0° direction 

45° flat sheet supported by a dimpled sheet oriented in the 45° direction 

ANALYSIS 

The assumptions used for the structural analysis of the multiwall sandwich are 
similar to those of reference 3. Since the structural behavior of flat sheets is 
well understood, the flat sheets were analyzed using standard techniques. Thus, 
the majority of the analytical effort in this paper is focused on the dimpled sheets 
because of their complex geometry. First, the extensional stiffnesses of the flat 
sheets in tension and compression are determined. The extensional stiffness of the 
dimpled sheet is separately determined. Then these extensional stiffnesses are used 
in simple beam theory calculations to estimate the bending stiffnesses of various 
configurations of multiwall sandwich. 
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Extensional Stiffness 

Flat sheets .- Because their primary function is to provide radiation barriers, 
the flat sheets, which are supported on the crests of the dimpled sheets, are made 
thin to reduce weight. Although they were assumed to be fully effective in tension, 
the thin flat sheets will buckle in compression at low stress. Appendix A shows 
the calculated compressive buckling stress of a flat sheet as a function of thick- 
ness. Since the buckling stress is very low for thin flat sheets, and thin sheets 
are likely to be prewrinkled, the flat sheets were assumed (as in ref. 3) to buckle 
immediately in compression and to have no significant postbuckling stiffness. 

Dimpled sheets.- In reference 3 the dimpled sheet was assumed to have the same 
inplane extensional stiffness as a flat sheet of the same thickness. To obtain a 
better estimate of the extensional stiffness, the dimpled sheet was analyzed using a 
finite element computer program, the SPAR structural analysis system (ref. 5). In 
this analysis the dimpled sheet was assumed to be linearly elastic and uniform in 
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thickness. The thickness used for the finite elements was obtained by multiplying 
the thickness of the dimpled sheet before forming t d by the ratio of the surface 
area of the flat sheet to that of the dimpled sheet. 

Two finite element models of the dimpled sheet were necessary because the geome- 
tries are different in the 0° orientation and the 45° orientation (see fig. 2) . 

Since the dimpled sheet has a repeating geometric pattern, only the small regions 
indicated in figure 2 were modelled. 

The finite element models, which represent the curved, three-dimensional surface 
of the dimpled sheet, are shown in figure 3. In both models, the tops of the dimples 
wete assumed to be parallel to the midplane of the dimpled sheet. Boundaries 1 and 2 
for the 0° direction model consist of two straight lines. One extends across the 
flat top of the dimple, and a second extends from the edge of the dimple top to the 
midplane. The diagonal connecting the intersection of boundaries 1 and 2 to the 
intersection of boundaries 3 and 4 consists of a straight line across the flat top of 

the dimpled sheet and a circular arc which is tangent to the midplane and connects to 

the edge of the dimple top. All four boundaries of the 45° model have the same 
shape as the diagonal of the 0° model. The diagonal connecting the dimples in the 
45° model has a shape similar to boundaries 1 and 2 of the 0° model. The perpendicu- 
lar diagonal is a straight line lying in the midplane. Triangular bending elements 
were used because the quadrilateral bending elements of SPAR require all four corner 
nodes to be coplanar for an accurate solution. To insure the symmetry of the struc- 
tural response, the triangular elements were superimposed as shown in figure 4, and 
each element was assigned only half of the material modulus so that the model would 
have the correct stiffness. As shown in figure 4, there is no node in the finite 
element model at the intersection of the diagonals of the overlapping triangular 
elements. The model oriented in the 0° direction, shown on the left-hand side of 

3, contains 100 nodes and 238 elements, and the model shown on the right-hand 

side of figure 3 and oriented in the 45° direction, contains 200 nodes and 
656 elements. 

For the model oriented in the 0° direction, boundaries 1 and 2 were treated as 
symmetry planes; i.e., translations and rotations about axes parallel to the plane of 
symmetry were constrained. Two different constraint cases were considered for bound- 
ary 3. In the first case, the boundary nodes were not allowed to translate in the 
x direction, and any lateral displacements were completely prevented. In the second 
case, a precalculated translation in the x direction was imposed on the boundary 
plane so that no net force existed normal to the boundary. The actual lateral con- 
straint provided to the dimpled sheet by the flat sheets in multiwall has not been 
determined, but it is bounded by these two cases. Displacement of the nodes on 
boundary 3 in the z direction was not allowed in either case. The model was loaded 
by a uniform displacement applied in the y direction to each node on boundary 4. 

For the model oriented in the 45° direction (fig. 3), the constraints applied 
to boundaries 1 and 2 were the same as those applied to boundaries 1 and 2 of the 
0° models; that is, they were treated as symmetry planes. The same two constraint 
conditions were applied to the x displacement of boundary 3 as for the 0° model to 
provide upper and lower bounds to the actual lateral constraint on the dimpled sheet 
in the multiwall sandwich. The nodes on boundaries 3 and 4 were not constrained from 
translation in the z direction but were constrained from rotation about the y 
and x axes, respectively. Also nodes on boundaries 3 and 4 were constrained from 
rotation about the z axis. The node at the corner of boundaries 1 and 4, the node 
®t the corner of boundaries 2 and 3, and the node midway between these two corners 
were fixed in the z direction. 
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The extensional stiffness of the dimpled sheet was determined from the applied 
displacements and the reaction forces calculated by using the SPAR finite element 
program. Average stress was calculated by dividing the sum of the reactions on the 
boundary opposite the applied displacement by the cross-sectional area t^L of the 
model. Th6 average strain was determined by dividing the applied displacement by the 
length of the model in the direction of displacement. The average stress was divided 
by the average strain to get an effective modulus E eff The extensional stiffness 
per unit width of the dimpled sheet is then E e f f t^ • Since the extensional stiffness 
of a flat sheet per unit width is Et £, the thickness of a flat sheet which will have 
the same extensional stiffness as the dimpled sheet is 


= Rt d 


where 


R 



eff 


E 


Thus R indicates the reduced extensional stiffness of a dimpled sheet. 

In addition to stiffness, the stress distribution in a dimpled sheet under an 
extensional loading was also determined by finite element analysis. The finite 
element models for both the 0° and 45° directions were refined by adding additional 
nodal points and triangular bending elements in areas of high stress to create the 
models shown in figure 5. The loading and boundary conditions remained the same. 
Stresses were calculated in the x and y directions and on the upper, middle, and 
lower surfaces at each node of boundaries 1 and 2 (fig. 5). These stresses were then 
divided by the overall average stress to obtain the stress concentration factors. 


Sandwich Bending Stiffness 

Once the extensional behavior of the flat and dimpled sheets is known, the 
moment of inertia and hence the bending stiffness of the multiwall sandwich can be 
easily calculated using standard engineering techniques. If the individual flat and 
dimpled sheets are assumed to have negligible bending stiffness compared with the 
total sandwich, the bending moment will be carried entirely by axial forces in the 
flat and dimpled sheets. By using some of the assumptions from reference 3, which 
are illustrated in figure 6, the moment of inertia of a multiwall sandwich can be 
determined by simple beam theory calculations. The sketch on the left of figure 6 
represents the actual multiwall sandwich loaded in bending. Ihe right-hand sketch 
illustrates the idealized structure in which flat sheets in compression have been 
removed entirely to simulate their buckling. Dimpled sheets have been replaced by 
flat sheets with their thicknesses reduced to represent the associated reduced exten- 
sional stiffness, and flat sheets in tension remain unchanged. Because the flat 
sheets in compression buckle and carry little load, the neutral bending axis is no 
longer at the geometric center of the structure. Calculations of moment of inertia 
for sandwiches of one, two, and three dimpled layers are shown in appendix B. Moment 
of inertia is used as a measure of bending stiffness because, for linear behavior, it 
is independent of material. 
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TESTS 


Tensile and bending tests were required to verify the assumption that thin, flat 
sheets in compression can be ignored and to determine experimentally the extensional 
stiffness of the dimpled sheet and the bending stiffness of various multiwall sand- 
wich configurations for comparison with analysis. The stiffness properties of 
interest, R for the dimpled sheet and I, moment of inertia, for the multiwall 
sandwich, are dependent only on geometry and should be independent of material. 
Therefore, any linearly elastic material can be used. Most of the tensile specimens 
and all the bending specimens were made of stainless steel because it was more 
readily available and more easily formed. However, since the present multiwall TPS 
tiles are made of titanium, some titanium tensile specimens were also tested. 


Sheet Tensile Tests 

Specimen fabrication. - Dimpled sheet specimens were fabricated from type 304 
stainless steel and 6A1-4V titanium, which were 0.004 and 0.003 in. thick, respec- 
tively. The dimpled sheets were formed from 12-in. by 24-in. flat sheets by using 
dies having evenly spaced but offset grids of protruding pins. The stainless steel 
sheets were cold formed; however, the titanium sheets were superplastically formed at 
an elevated temperature. Dogbone-shaped, stainless steel specimens, 10 in. long and 
1 .5 in. wide at the test section, and 3-in. by 1 2-in. rectangular titanium specimens 
were cut from the dimpled sheets in both 0° and 45° orientations with an electronic 
discharge machine (EDM). 

Test setup and procedure . - The specimens were loaded in a universal test 
machine, with the machine crosshead moving at a constant rate of 0.1 in/min. Load 
was measured by a simple reading from a load cell in the machine; however, diffi- 
culties were encountered in attempting to measure strain. Standard strain gauges 
could not be used to measure average strain because of the uneven surface of the 
specimens. This uneven surface also made it difficult to properly attach available 
extensiometers to the specimen. Because the specimens were tested at low loads and 
the test machine was designed for much higher loads, the strain was assumed to be 
the displacement of the machine crosshead divided by the original length of the 
specimen between the grips. Thus, input from the load cell and from the motion of 
the crosshead were used to record load-deflection curves on the chart recorder of 
the test machine. 


Sandwich Bending Tests 

Specimen fabrication .- Bending test specimens, shown in figure 7, were fabri- 
cated from 12-in. by 24-in. sheets of type 304 stainless steel brazed together. The 
flat sheets were 0.002 in. thick, and the dimpled sheets were 0.004 in. thick before 
dimpling. The dimpled and flat sheets were assembled to make multiwall panels of 
one, two, and three dimpled layers. Holes were drilled, and pins were inserted 
through the corners of the panels to assure alignment of the dimpled sheets. The 
assembled dimpled and flat sheets were brazed at 1470°F for 2 to 4 min in a vacuum 
furnace. For each of the one-, two-, and three-layer configurations, specimens 
(2 in. by 11 in.) were fabricated with dimpled sheets oriented in the 0° direction 
and in the 45° direction. 

Test setup and procedure . - A sketch of a standard four-point bending test 
arrangement is shown in figure 8. The specimen was supported by 1/2-in-diameter 
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steel rods spaced 10 in. apart in a metal framework, which rested on top of a 
load cell. The load was applied through two 1 /2-in-diameter steel rods, spaced 
6 in. apart and attached to the moving machine head of a universal test machine. 

The deflection was measured from the center of the specimen relative to the moving 
machine head by using a direct current displacement transducer (DCDT) • Load- 
deflection curves were recorded by a flat-bed plotter from the output of the DCDT 
and the load cell. 

Multiwall specimens of one, two, and three dimpled layers with dimpled sheets 
oriented in both the 0° and 45° directions were tested in four-point bending tests. 
Each specimen was loaded to a point in the linear range of the load-deflection curve 
and unloaded two times prior to loading to failure. Thus, at least three load- 
deflection curves were recorded for each specimen. The specimens were first tested 
intact to obtain their experimental bending stiffnesses. Then, for each of the one- 
and two-dimpled-layer specimens, the outer flat sheet which had been in compression 
was cut at each dimple pitch so that it could carry no load. These specimens were 
retested to determine how much the flat sheets in compression contributed to flexural 
stiffness of the multiwall sandwich. Then the other outer flat sheet of each two- 
dimpled-layer specimen was also cut so that only the two dimpled layers remained to 
carry the bending load. These two-dimpled-layer specimens were retested, and the 
data were used to determine the extensional stiffness of the dimpled sheet by the 
method in appendix C. 


RESULTS 

Extensional Stiffness 

Flat sheets. - One of the key assumptions in the multiwall design approach pre- 
sented herein and in reference 3 is that flat sheets in compression buckle at a low 
load and have a negligible contribution to flexural stiffness. Results of bending 
tests on one- and two-dimpled-layer specimens indicate this assumption is valid when 
the flat sheets are 0.002 in. thick. Figure 9 shows how the flat sheets in compres- 
sion affect the moment of inertia of multiwall bending specimens. The moment of 

inertia, which is proportional to bending stiffness, is shown for four different 
multiwall configurations. The method used to calculate the experimentally determined 
moments of inertia from the load-deflection curves is presented in appendix D. The 
unshaded bars in figure 9 represent the moment of inertia of intact specimens, and 

the shaded bars represent the moment of inertia of specimens with the flat sheets in 

compression effectively removed. The moments of inertia are nearly the same for the 
intact specimens and the cut specimens. Thus, it appears to be a valid assumption to 
ignore the 0.002-in-thick flat sheets in compression when calculating bending stiff- 
ness of multiwall. 

If, however, these flat sheets are thicker than 0.002 in., they may have a 
significant contribution to the bending stiffness. The calculations described in 
appendix E and summarized in figure 10 indicate the effect that thicker flat sheets 
would have on the behavior of a multiwall sandwich in bending. Figure 10 shows the 
predicted deflection curves which would be obtained in four-point bending tests of 
two-dimpled-layer specimens with flat sheets of different thickness. These curves 
were calculated by assuming that the flat sheets buckle at stresses predicted in 
appendix A and that the flat sheets have no post buckling stiffness. The initial 
slope of each curve reflects the flexural stiffness before the flat sheets buckle, 
and the final slope reflects the bending stiffness after the flat sheets buckle. The 
horizontal portion of each curve represents the load at which the flat sheets reach 
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their compressive buckling stress. The dashed line is a load-deflection curve 
measured for a specimen with 0.002-in. flat sheets. The slope of the measured curve 
is close to that predicted for the specimen with buckled flat sheets at low loads. 
These calculations indicate that the flexural behavior of the multiwall sandwich is 
sensitive to the flat-sheet thickness and that the assumption that flat sheets buckle 

and carry no further load in compression is not valid as flat-sheet thickness is 
increased. 

A comparison of the approach taken in this paper with that of reference 6 illus- 
trates the importance of the flat-sheet behavior in compression. The multiwall bend- 
ing specimens tested in reference 6 had 0.003-in-thick flat sheets. The flat sheets 
remained unbuckled for a measurable portion of the load-deflection curve, and bending 
stiffnesses in reference 6 are calculated from the initial slope of the load- 
deflection curve, where the flat sheets were not buckled. As a result, the stiff- 
nesses presented in reference 6 are much higher than would be predicted by the 
analysis presented in this paper. 

Dimpl ed sheets . - Another key assumption in the multiwall design approach is 
that the structure exhibits linear behavior. Since high stresses can cause material 
nonlinearity, stress concentration factors were calculated for the dimpled sheet. 
Figure 11(a) shows stress distributions for the 0° orientation of the dimpled sheet 
calculated with the SPAR structural analysis system. Stresses are normalized by the 
average stress in the direction of loading, so that the curves shown represent stress 
concentration factors. Stress concentrations are shown on the boundaries, which 
are the locations of maximum stress on the model. Curves are shown for the upper, 
center, and lower surfaces of the bending elements. The large differences between 
stresses on the upper and lower surfaces indicate bending. The locations and signs 
of the bending stresses correspond to observed deformation of tensile test specimens. 
The large bending stresses along the x axis correspond to the observed flattening 
of the dimple. Along the x = L boundary, the stress concentrations indicate there 
is little bending. However, the maximum principal stresses, which are within 10° of 
the y-direction stresses, indicate that much of the load is being carried along this 
path. The stress concentrations are high along both boundaries and exceed 4 at 
several locations. 

Similarly, stress concentrations were calculated for both constraint cases of 
the dimpled sheet oriented in the 45° direction. Figure 11(b) shows the stress con- 
centrations for the case of an allowed lateral displacement. This constraint case 
approximates the behavior of the dimpled sheet tensile test. As in figure 11(a), 
curves are shown for the upper, center, and lower surfaces of the bending elements. 
These curves indicate that there are very large stress concentrations due to bending 
along both boundaries of the model. The stress concentrations due to bending exceed 
25 at the lower right-hand corner of the model. The reason for these large stress 
concentrations is apparent. The dimpled sheet tends to stretch in the direction of 
loading and to narrow or contract transverse to the direction of loading. Both of 
these deformations are resisted primarily by localized bending, which results in 
large bending stresses. 

Stress concentrations for the case of no lateral displacement of the 45° dimpled 
sheet are shown in figure 11(c). Although the stress concentrations have a distribu- 
tion similar to that in the previous example, the magnitudes of the stress concentra- 
tions are less than one quarter those of the 45° dimpled sheet with lateral displace- 
ment. The maximum stress concentration factor is about 6 and, again, is located in 
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the lower right-hand corner of the model. This dramatic decrease in bending stress 
illustrates how the structural behavior of the dimpled sheet is strongly dependent on 
the lateral constraint. 

Stress-strain curves measured from tensile tests of both stainless steel and 
titanium dimpled sheets are shown in figure 12. The stress shown is the average 
stress obtained by dividing the load by the cross-sectional area. The average strain 
was determined by dividing the machine-head motion by the original length of the 
specimen. For a flat sheet tested in the same machine over the same load range as 
the dimpled sheet, the strain was measured by a strain gauge and found to be approxi- 
mately half that determined from machine-head motion. Therefore, the strain should 
be approximately half of that shown in figure 12. Although the strain is not accu- 
rate, the curves in figure 12 illustrate the characteristic behavior of the dimpled 
sheets and the differences in behavior between sheets oriented in the 0° and 
45° directions. 

Both the stainless steel and titanium dimpled-sheet specimens show similar 
behavior in the 0° direction. The specimens behave linearly until the material 
yields and the stiffness is reduced. The dimples tended to flatten, and there was 
no noticeable narrowing of the specimen. The predicted average stresses at which 
local material yielding occurs for the 0° dimpled sheet are indicated on figure 12. 
This stress was calculated by dividing the material yield stress (stainless steel, 

42 000 psi; titanium, 130 000 psi) by the largest stress concentration factor (4.9) 
from figure 11(a). These predicted stresses are close to the stresses where the 
dimpled sheets begin to show nonlinear behavior. The stress-strain curve for the 
0° direction indicates that the linear finite element analysis will be adequate to 
predict the behavior of the 0° dimpled sheet at low loads. 

The dimpled sheet exhibits completely different behavior in the 45° direction. 
The stress-strain curves of the 45° dimpled sheet (see fig. 12) show that at low 
loads, the dimpled sheet has little stiffness, but as the load increases, the stiff- 
ness also increases. This increased stiffness is caused by the specimen changing 
shape under load. Figure 13 shows a stainless steel, 45° dimpled-sheet specimen 
before and after deformation under a tensile load. Each vertical line of dimples is 
pulled into a corrugation, and the specimen is narrowed significantly, as indicated 
by the dashed lines. Although the dimpled-sheet specimen is obviously deformed con- 
siderably beyond the range of linear behavior, the figure illustrates the type of 
deformation which occurs. 

To understand this observed nonlinear behavior, it is helpful to consider a 
dimpled sheet loaded in tension as a series of parallel beams loaded in tension. In 
the 0° direction, some of these beams are straight and provide a straight-through 
load path, while others are wavy. The straight beams act as rods and carry the load 
in tension, but the wavy beams carry the load primarily through bending. The 
straight beams are therefore much stiffer and, consequently, carry more of the load. 
As the load increases, the strain increases, and the curved beams begin to straighten 
out, which results in flattening of the dimples. 

For the 45° direction, it is helpful to consider the dimpled sheet as a lattice 
of straight and wavy beams crisscrossed at ±45°. There is no straight-through load 
path. That is why the dimpled sheet has less stiffness in the 45° direction than in 
the 0° direction. When this lattice is loaded in tension, the straight beams tend to 
rotate to align themselves with the load. This induces lateral contraction, which 
causes the dimpled sheet to narrow significantly and to fold into corrugations with 
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crests running parallel to the direction of the applied load. The resulting corru- 
gated shape significantly reduces the curvature of the wavy beams and correspondingly 
increases the overall stiffness of the dimpled sheet. 

The importance of lateral boundary constraint becomes apparent when considering 
the 45° dimpled sheet as a lattice of rods. When the lateral edge is unconstrained, 
as in the tensile test, the lateral contraction is resisted only by localized bend- 
ing of the dimpled sheet; however, when the edge is constrained, the lateral contrac- 
tion is resisted at the boundary, and the localized bending is inhibited. Since 
this shape change is reduced, the behavior of the dimpled sheet should be more nearly 
linear. In the multiwall sandwich, a partial lateral constraint is imposed on the 
dimpled sheet by the flat sheets and other dimpled sheets. Consequently, the 
45° dimpled sheet may show more nearly linear behavior in the multiwall sandwich than 
in the tensile tests. 


Although the stress-strain curves of both the stainless steel and titanium 45° 
dimpled sheets (see fig. 12) show the same trends, their shapes are somewhat differ- 
ent. The two most likely causes are different thicknesses and different materials. 
The stainless steel dimpled sheet was 0.004 in. thick and the titanium sheet was 
0.003 in. thick; however, analysis indicates that this difference in thickness would 
not cause the behavior to differ as much as observed. Stainless steel has a yield 
point which is one-third that of titanium. Considering the high initial stress con- 
centration factors predicted by the finite element analysis, some localized yielding 
is likely to occur in the stainless steel before it would occur in the titanium. 

This localized yielding may enable the stainless steel to deform more quickly into 
the stiffer shape. The offsetting factors of increased stiffness due to shape change 
and decreased stiffness due to localized yielding may result in the shape of the 
stainless steel stress-strain curve* 

Because the tensile test specimens did not have lateral boundary conditions 
representative of the multiwall sandwich, bending tests of sandwich specimens were 
used to measure the extensional stiffness of the dimpled sheet. The method used to 
calculate the extensional stiffness of the dimpled sheet from the measured load- 
deflection curves is presented in appendix C. In figure 14 the dimpled sheet exten- 
sional stiffnesses obtained from bending tests are compared with those predicted 
using the SPAR finite element program. The extensional stiffness is indicated in 
terms of R, which is the ratio of the thickness of a flat sheet t f to that of a 
dimpled sheet t d with the same extensional stiffness. Therefore, by definition, 
for a flat sheet, R = 1 . The open circles represent the extensional stiffness of 
the 0° dimpled sheet, and the triangles represent that of the 45° dimpled sheet 
obtained from bending tests. Each symbol represents a single specimen. The curves 
on the left-hand side of the figure represent the extensional stiffnesses of the 
dimpled sheet calculated for various dimpled-sheet thicknesses t, and a constant 
dimple height h^. The curves on the right-hand side represent extensional stiff- 
nesses calculated for various dimple heights h d and a constant dimpled-sheet thick- 
ness t d . The upper narrow band represents the calculated extensional stiffness for 
the 0° direction; the lower wide band represents the calculated extensional stiffness 
for the 45° direction. The experimental and predicted results agree well in that 
they both show the extensional stiffness of dimpled sheet in the 0° direction to be 
about 30 percent that of a flat sheet, and in the 45° direction to be about 10 per- 
cent that of a flat sheet. For the 0° direction the predicted stiffness falls within 
the experimental scatter, and for the 45° direction the test data fall within the 
predicted range. 


10 



Bending Stiffness 


In figure 15 moments of inertia determined from the bending tests of multiwall 
specimens with one, two, and three dimpled layers are compared with moments of 
inertia predicted using values of R calculated by the finite element analysis. The 
shaded bars represent experimental values, and the unshaded bars represent analyti- 
cal values. The experimental values were obtained from the measured load-deflection 
curves by the method given in appendix D, and the predicted values were obtained 
by the method given in appendix B. Since a range of R was predicted by the finite 
element analysis because of lateral displacement boundary conditions, there is a 
corresponding range of calculated values of moment of inertia indicated by the 
arrows. The experimental values fall well within this range. Considering the 
complex behavior of the dimpled sheet observed during tensile tests, the agreement 
between analysis and experiment is surprisingly good. Although the dimpled sheet 
behaves nonlinearly in tension, multiwall sandwich specimens with dimpled sheets 
oriented in the 0° or 45° direction exhibited linear behavior in bending. These 
results indicate that the bending stiffness of the multiwall sandwich can be deter- 
mined by a relatively simple engineering approach. However, additional test data 
are required to confirm this conclusion for other multiwall geometries and fabrica- 
tion methods. 


CONCLUDING REMARKS 

An analytical and experimental study was carried out to understand the struc- 
tural behavior of multiwall sandwich. The two primary components of multiwall sand- 
wich, flat sheets and dimpled sheets, were analyzed separately to determine behavior 
in tension and compression. Once this behavior was known, the moment of inertia of 
the multiwall sandwich, which is proportional to bending stiffness, was determined 
using standard engineering calculations. Test data were obtained to evaluate the 
accuracy of these bending stiffness calculations and to evaluate the assumptions used 
in the analysis. 

Comparison of the test data with the calculated values verifies the assumption 
that the flat sheets which are loaded in compression during bending will buckle, 
carry little load, and can be ignored. However, the flat sheets in these specimens 
were only 0.002 in. thick, and calculations indicate that thicker flat sheets will 
carry a more significant load before buckling. Consequently, the assumption that 
flat sheets in compression can be ignored becomes less valid if the thickness of the 
sheets is increased. 

Both the finite element analysis and the test data show that the dimpled sheet 
has significantly less stiffness than a flat sheet of the same thickness. Also, 
analysis and experiment show that the dimpled sheet has directional properties and 
is significantly stiff er in the 0° direction than in the 45° direction. Tensile 
tests show that in the 0° direction, the dimpled sheet behaves linearly until mate- 
rial yielding occurs; however, in the 45° direction, the dimpled sheet exhibits 
nonlinear behavior as a result of its significant shape changes. Finite element 
analysis shows that high local bending stresses occur in the tensile specimens and 
that the associated shape changes are strongly influenced by the lateral-edge 
constraints. The dimpled sheets in the multiwall sandwich are not free from lateral- 
edge constraint. Therefore, the behavior of the 45° dimpled sheet in the multiwall 
sandwich may be much more nearly linear than that observed in the tensile tests. Two 
constraint cases were used in the analysis to bracket the actual lateral constraint 
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imposed on the dimpled sheet in the multiwall sandwich. Results of this analytical 
approach compare well with test data. Both show that the extensional stiffness of 
the dimpled sheet in the 0° direction is about 30 percent of that for a flat sheet, 
and that in the 45° direction, it is about 10 percent. 

The predicted multiwall— sandwich bending stiffness agrees surprisingly well 
with experiment for all sandwich configurations tested. Although the dimpled sheet 
behaves nonlinearly in tension, multiwall-sandwich specimens with dimpled sheets 
oriented in the 0° or 45° direction exhibited linear behavior in bending. The close 
agreement between analysis and experiment indicates that the method of calculating 
bending stiffness presented in this paper is valid for this linear range of behavior. 
However, for other multiwall geometries and fabrication methods, this method may be 
less accurate. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 11, 1983 



APPENDIX A 


BUCKLING OF FLAT SHEETS IN COMPRESSION 

The flat sheets which are in compression during bending are assumed to buckle at 
a low stress and, consequently, to carry no load. The calculations in this appendix 
estimate the buckling stress of the flat sheet. Since only uniaxial bending is con- 
sidered, buckling of the flat sheets is approximated using column theory. The buck- 
ling load should be lower than predicted because the boundary condition is actually 
between simply supported and clamped, and the foil-gauge sheet will likely not be 
perfectly flat. 

The following sketch represents a flat sheet in uniaxial compression. Each dot 
represents a point at which the flat sheet is supported by a dimpled sheet oriented 
in the 0° direction. The curved lines represent alternate crests and troughs, which 
are the buckling shapes observed during the bending tests. The area enclosed by 
dashed lines was modelled as a clamped column, buckling in the second mode. 


0° Orientation 



From column buckling theory, 


(P ) 
cr 0 


8.1 8ti 2 EI 
(\f2s) 2 


where (P cr ) 0 o 


( °cr ) 0° bt f ' and 1 



(A1 ) 
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By substitution, 


( 

cr 


2 2 
8.1 8u Et 


24s ‘ 


(A2) 


A similar sketch is shown below for a flat sheet supported by a dimpled sheet 
oriented in the 45° direction. In this case, the observed buckling shape was much 
simpler. The flat sheet buckles upward between each lateral row of dimple peaks. 
Consequently, the area enclosed by the dashed lines was treated as a clamped column 
buckling in the first mode. 


45° Orientation 



Section B-B 


From column buckling theory. 


{p 

cr 45 


4 m 2 El 

2 

s 


where (P ) „_ 0 
cr 45 


(a cr ) 45 obt f' and 1 



(A3) 


By substitution, 


( a )._o = 

cr 45 


m 2 Et f 2 


( A4 ) 
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The following figure was generated from the two expressions for a cr developed in 
this appendix. The figure shows that the buckling stress is very low for the 
0.002-in. flat sheets which were tested, but that the buckling stress increases 
quadratically as flat-sheet thickness increases. 


Supported by 45° dimpled sheet 
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APPENDIX B 


CALCULATION OF MOMENT OF INERTIA OF MULTIWALL BENDING SPECIMENS 

This appendix contains calculations which predict the moment of inertia of the 
bending test specimens of one, two, and three dimpled layers. The flat sheets 
loaded in compression are assumed to buckle and carry no load. The dimpled sheet 
is modelled as a flat sheet with a reduced thickness lying along the midplane of the 
dimpled sheet. The reduced thickness represents the reduced stiffness associated 
with the dimpled sheet. The moment of inertia of the resulting model is then cal- 
culated using the parallel axis theorem. 


One-Dimpled-Layer Multiwall Sandwich 

The idealization of a one— dimpled— layer multiwall sandwich loaded in bending is 
shown in the following sketch. 


Unit width cross section 



Modelled as 




The neutral bending axis y is calculated by the following equation: 


y = JL5SZ 

3 E EA 


For the one-dimpled-layer multiwall. 


(B1 ) 


y = 


(0)t f + Rt dl~ 

*f + Rt d 


( B2 ) 


Once y is known, the moment of inertia I can be calculated from 


I = T (I. +d.A.) 

1 11 

i=1 


(B3 ) 
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where n is the number of sheets , and d is the distance from the neutral bending 
axis. The moment of inertia of each flat or dimpled sheet 1^ is assumed to be 
negligible compared with the moment of inertia of the total sandwich. 


I. « I 
1 


Thus, 


I = t f (y ) 2 




( B4 ) 


Two-Dimpled-Layer Multiwall Sandwich 

The calculations are similar for a two-dimpled-layer sandwich. The following 
sketch shows the idealized structure. 



Again, the location of the neutral bending axis is calculated from 
equation (B1 ) . 


y 


(0)t f 


+ Rt d\2 

+ 2Rt 
f d 



( B5 ) 


where R is from figure 14. From equation (B3), the moment of inertia I is 


I = t f (y) 


2 


+ Rt 

d 




2-i 


(B6) 
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Three-Dimpled-Layer Multiwall Sandwich 

Similarily, for a three-dimpled-layer multiwall sandwich, the idealized struc- 
ture is shown in the following sketch. 




A = Rt 


d 



The neutral bending axis is 


y = 


h 3h 5h 

t f (0 + V + R t d l + 

2t + 3Rt _ 
f d 


Simplifying yields 


y 


h d( fc f + 2 Rt d) 

2t c + 3Rt , 
r d 


( B7 ) 


If y < h d ' then th e second-from-bottom flat sheet is loaded in compression and, 
thus, is effectively removed, and y must be recalculated. 

If y < hj, then 


y* 


V 0) + 2 Rt d h d 

t. + 3Rt . 
t d 


(B8) 
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and 


I 


t f (y *) 2 + Rt d 





(B9) 


where y* is the neutral bending axis calculated by assuming that all flat sheets 
except the bottom sheet have buckled. 



(BIO) 
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APPENDIX C 


CALCULATION OF R FROM EXPERIMENTAL DATA 

This appendix contains the method of calculating R, the ratio of the stiffness 
of a dimpled sheet to that of a flat sheet with the same thickness, for a dimpled 
sheet using data from a four-point bending test. The moment of inertia of the two- 
dimpled-layer specimen, with upper and lower flat sheets cut to totally remove their 
contribution, was calculated from the load-deflection curve by the method given in 
appendix D. If the geometry of the dimpled sheet (sketch below) and the moment of 
inertia of the specimen are known, R for the dimpled sheet can be calculated in the 
following manner. 



The neutral bending axis y is 



The moment of inertia I is 



Thus 


R 


21 


fc d h d 


(Cl ) 


(C2 ) 


(Use the experimentally determined value of I to obtain R.) 
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CALCULATION OF MOMENT OF INERTIA FROM FOUR-POINT BENDING TESTS 

This appendix contains the method used to calculate the experimental moment of 
inertia from the measured loads and deflections of the test specimen. Deflections 
were measured at the center of the specimen, with respect to the points where the 
load was applied. The inner section of the specimen, between the points where the 
load was applied, was analyzed as a separate beam loaded in pure bending. An equa- 
tion was developed relating the measured slope of the load-deflection curve to the 
moment of inertia of th6 specimen. 


P 

2 

4 — a— i 




V I I ) M o 

V I -X 



The central section of the specimen shown above can be treated as a separate 

. ( ap\ 

beam of length 1 loaded by pure bending IM Q = — . 

The moment-curvature relation for a beam is 


.2 

d w 

El = M 

dx 


M o 


aP 

2 


where w = Transverse displacement. Integration gives 


„ dw aP , _ 
El — = — X + C 
dx 2 1 


Applying the boundary condition to solve for yields 


dw 4 = 0 = + c 

dx V2 / 4 1 


C 


1 


aPA 

4 


( D 1 ) 
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Integrating again gives 


.. aPx aPl „ 

EIw = — — x + C 

4 4 2 


Applying the boundary condition to solve for C 2 yields 


"0 ■ 0 - C 2 


w = 


aPx aPlx 
4EI “ 4EI 


( D2 ) 


At x = 1/2, 6 is the measured deflection, 


w 


= -6 = 


aPl 


a PI 


1 6EI 8EI 


6 = 


a PI 
16EI 


( D3 ) 


Solving for I gives 


I = 


a l 2 _P 
16E 6 


( D4 ) 


where P/6 is the measured slope of the load-deflection curve, 
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EFFECT OF FLAT-SHEET THICKNESS ON PREDICTED LOAD-DEFLECTION CURVES FOR 

FOUR-POINT BENDING TESTS 

The calculations in this appendix show how the use of thicker flat sheets may 
affect the flexural behavior of multiwall sandwich. As an example, load-deflection 
curves were predicted for a four -point bending test of a two-dimpled-layer multiwall 
sandwich with dimpled sheets oriented in the 0° direction. The shape of the curve 
was predicted by calculating the slope of the curve before the flat sheet buckled, 
the load at which buckling occurred, and the slope after the flat sheet buckled. 

To calculate the slope before buckling, the moment of inertia must be calcu- 
lated. The following sketch shows an idealized two-dimpled-layer sandwich loaded in 
bending with no flat sheets buckled. 



From equation (B3), 


I,, = z (I. + d. 2 a. ) 

bb “ l 11 

i=1 


where n = Number of sheets and I bb = I before buckling of flat sheets. Assuming 
I ± « If 



To obtain the slope before buckling, equation (D4) is used. 

= 2 
J hh 1 6E 6 
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Thus, the slope of the load-deflection curve before flat sheets buckle is given by 


M 16EI bb 
' 6 'bb aZ 2 


(E2) 


To obtain the load at which buckling occurs , equation (A2) for buckling stress is 
used. For a flat sheet supported by a 0° dimpled sheet. 


(a 

cr 


) 


0 ° 


2 2 

8 .1 Qn Etf 
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And from simple beam theory. 


O 


M(y - y) 

^b 



( E3 ) 


For 0 


(cx 


cr 


0 °' 



2 I bb (g crV 

ah , 


= Load at which flat sheet buckles 


(E4) 


To determine the slope of the curve after buckling, the moment of inertia must be 
calculated from equation (B5). 


2Rt,h 

- a d 

y “ t. + 2Rt , 

t d 

From equation (B6), 


I 


after buckling = 


ab 


= t f (y ) 2 + 




2^ 


From equation (D4), 


p\ ,6Ei ,b 

T ) = 5 — = slope of curve after buckling of flat sheets 

/ab af 

These equations were used to generate the curves shown in figure 10. 


(E5) 
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Figure 1 Titanium multiwall sandwich 


Dimpled sheet 



HI Region modelled for 0° 
“ orientation (1/4 dimple) 



Region modelled for 45° 
orientation (1/2 dimple) 


O Node above midplane 


• Node below midplane 



Figure 2.- Regions of dimpled sheet modelled for structural analysis. 
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0° direction 



Figure 3,- SPAR finite element models of dimpled sheet. 


Finite element model 



Figure 4.- Illustration of superimposed triangular elements. 
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0° orientation 45° orientation 


Figure 5.- Finite element models used to calculate stress concentration factors. 


Structure 


Idealized structure 




Dimpled sheets in compression 



Neutral 

bending 

axis 



Flat sheet in tension 


Figure 6.- Illustration of simplifying assumptions used to calculate bending stiffness of multiwall sandwich. 





Figure 8.- Four-point bending test setup. 
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Figure 9.- Effect of flat sheets in compression on moment of inertia of multiwall 

sandwich, t^ = 0.002 in. 




Figure 10.- Predicted load-deflection curves for four-point bending tests on multiwall sandwich. 
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(a) 0° direction. 

Figure 11.- Stress concentration factors in dimpled sheet. 
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(b) 45° direction; lateral displacement allowed. 


Figure 11.- Continued 
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(c) 45° direction; lateral displacement constrained. 


Figure 11.- Concluded. 
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Figure 12.- Stress-strain curves of dimpled sheets (determined from 

machine-head motion) . 
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Figure 15.- Effective bending moments of 
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